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ABSTRACT 

We have used the Submillimeter Array for the first interferometric CO imaging toward the starburst-Seyfert 
nucleus of the southern barred spiral galaxy NGC 1365, which is one of the four galaxies within 30 Mpc that 
have Lg-iooo^m > 10 n L Q . Our mosaic maps of 12 CO, 13 CO, and C ls O(J=2-l) emission at up to 2" (200 pc) 
resolutions have revealed a circumnuclear gas ring and several CO clumps in the central 3 kpc. The molecular 
ring shows morphological and kinematical signs of bar-driven gas dynamics, and the region as a whole is 
found to follow the star formation laws of Kennicutt. We have found that some of the gas clumps and peaks in 
CO brightness temperature, which we collectively call CO hotspots, coincide with the radio and mid-infrared 
sources previously identified as dust-enshrouded super star clusters. This hotspot-cluster association suggests 
that either the formation of the most massive clusters took place in large molecular gas concentrations (of 
S mo i ~ 10 3 Mq pc~ 2 in 200 pc scales) or the clusters have heated their ambient gas to cause or enhance the 
CO hotspots. The active nucleus is in the region of weak CO emission and is not associated with distinctive 
molecular gas properties. 

Subject headings: galaxies: starburst — galaxies: ISM — galaxies: star clusters — galaxies: active — galax- 
ies: individual (NGC 1365) 



1. INTRODUCTION 

Central regions of spiral galaxies sometimes host vigorous 
star formation called (circum)nuclear starbursts. This phe- 
nomenon forms stars at a higher rate and more preferentially 
in com pact massive clusters than a non-s tarburst environment 
(e.g.,|jelesco et alll993tbarth et all 19951) . Since stars form 
from molecular gas, properties of the molecular gas are ex- 
pected to control the starbursts. The study of such properties 
in nearby starbursts tells us about physical processes involved 
in the phenomenon, and the information serves as a basis for 
modeling starbursts throughout galaxy evolution. 

The typical extent of a (circum)nuclear starburst in our 
vicinity is < 1 kpc, which is only 20" for a distance of 10 
Mpc. High angular resolution is therefore essential for star- 
burst studies. In particular, high resolution imaging of molec- 
ular gas with millimeter interferometers has been a power- 
ful tool in the field (e.g., iTo gee et aTj |2005). Until very re- 
cently, however, such arrays had been operated only at the 
mid-latitudes in the northern hemisphere, leaving many south- 
ern galaxies unattainable. 

We have now observed NGC 1365, a starburst-Seyfert com- 
posite galaxy in the southern sky, using the Submillimeter Ar- 
ray (SMA) 6 . This is part of our SMA survey of southern star- 
bursts, aiming to explore molecular gas and dust in some of 
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the most important but undeservedly less observed starburst 
galaxies in the southern sky. NGC 1365 is one of our targets 
because it has the circumnuclear starburst typical of barred 
spirals, and because the infrared luminosity of the starburst 
is among the largest in non-merging galaxies, bridging more 
luminous starbursts in distant mergers and less luminous ones 
in undisturbed spirals. Other galaxies reporte d so far from our 
survey are M8 3, NGC 253, and NGC 3256 dSakamoto et alJ 
l2004ll2006aibT) . 

We have observed CO emission to trace the molecular gas. 
Our resolution of 2" is a significant improvement over the best 
previous ones of 5"-10". In addition, the three isotopic lines 
that were simultaneously imaged, i.e., the J=2-l transitions of 
12 CO, I3 CO, and C I8 0, provide information about the physi- 
cal properties of the molecular gas. Among other things, we 
have discovered bright CO spots, or CO hotspots, that are as- 
sociated with young massive star clusters in the starburst. We 
discuss their properties and relation to the starburst. 

The goals of this paper are (1) to report the detailed struc- 
ture, kinematics, and properties of molecular gas in the central 
region of NGC 1365 on the basis of the new high-resolution 
data, (2) to search for the relation between circumnuclear star- 
burst and the gas properties, and (3) to see if there is anything 
special in the molecular gas in the vicinity of the Seyfert nu- 
cleus. The next section further introduces NGC 1365, and ^3] 
describes our SMA observations and data reduction. The ob- 
servational results are presented in 5j4] and are discussed in §5\ 
- §Q A summary is given in ^8] 

2. NGC 1365 

NGC 1365 is an archetype of barred spiral galaxies and 
hosts a circumnuclear starburst and a Seyfert nucl eus. A com - 
prehensive review of the galaxy is given in Lindblad ( 1999). 
A short summary of the properties of this galaxy is provided 
below and in Tabled We a dopt a Cepheid-based distance of 
17.95 Mpc (iFreedman et al.ll200ll) . At this distance, 1" cor- 
responds to 86 pc and 1 kpc to 11 "5. The galaxy has a long 
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stella r bar, spanning about 3' on the sky (Regan & Eh negree"r] 
119971) or 1 8 kpc on the galaxy plane. The bar has a pair of dust 
lanes in its front edge, assuming that the spiral arms are trail- 
ing in the galaxy. The assumption makes the northwest side 
of the galaxy near to us. The dust lanes curve inward near 
the galactic center to suggest an oval ring around the nucleus 
(see, e.g., the optical images in NED). 

Active star formation in the circumnuclear region was 
first noticed as optical hotspots by Morgan ( 1958J and 
ISersic & Pastorizal {l967). Later optical imaging using the 
Hubble Space Telescope (HST) reveal ed dozens of brigh t 
compact clusters around the dust ring ( Kristen et al. 1997). 
The brightest of them, after extinction correction, is 300 
times brighter than the most luminous glob ular cluster in 
our Ga laxy. Mid-IR (10-13 /jm) imaging by Galliano et al. 
(2005) also revealed several compact sources, presumably 
compact clusters, in the region. These clusters have been 
identified as the so-called super star clusters, or young clus- 
ters as massive as ^10 6 M Q containing thousands of OB stars. 
The infrared luminosity of NGC 1365 is Ls-iooo^m = 10 11 00 
Lq dSanders et all2003l) . which makes it one of only four ob- 
jects having > 10 11 L Q within 30 Mpc 7 . At 15 /im, 71% of the 
total fl ux of the galaxy originates from within the central 43" 
( iRoussel et alJl200"ll) . Thus most of the infrared luminosity 
probably comes from inside the starburst region of less than 2 
kpc radius. 

In the radio, where dust obscuration is minimal, a circum- 
nuclear ring of about 20" x 9" containing several co mpact 
(< 1" ) peaks have been found at 2-20 cm (Sandavis^j^D 
M9821 119951 ISaikia et all Il994[ iForbes & Norrisl 119981 
Morg anti et aflll999t iThean et al.ll2000t) . Some of the radio 
and the mid-IR sources coincide. The brightest radio hotspots 
have a power of P6cm~10 20 W Hz" 1 for each. In other words, 
they are more than 100 times as luminous as Cas A, are 
comparable to the brighte st radio sources in M82 and NGC 
253 ilKronberg et1dll985lluivestad & Antonuccill9'97l) . and 
are also co mparable in radio power with radio supernovae 
dWeiler et alJ2002l) . The center of the galaxy is als o bright in 
molecular lines ("e.g.. lSandavist etaflll995L 11999). Previous 
observatio ns found bar-driven twin peaks of mol ecular gas in 
the region (Sandavist et al. 1999; Ott et al. 2005), and higher 
resolution observations fully resolving the gas structure have 
been needed. 

The Seyfert 1.5 nucleus of the galaxy is the driver, or 
at least the ionizing source, of a biconical plasma o utflow 
dPhillips et al.ll983(lKristen et alll997UVeilleux et alJ2003l) . 
but the AGN is of low-luminosity and a s mall con tribu- 
tor to the luminosity of the ce ntral region (llyomoto et"aTI 
H9971 iKomossa & S chulzl fT998t LStevens et al.lll999l) . The 
mid- and far-IR spectral indexes of the galaxy, a(60,25) = 
-2.2 and a(100,60) = - 1.1 fo r S x oc X~ a and the IRAS 
fluxes in Sanders et a D d2003), 

are not those of 'warm' 
AGN-dominate d sources but those o f H II-type galaxies 
dde Griip et al.lll985t iMilev et al.l [l985). The position of the 
nucleus in Tablets the average of optical, near-IR, and radio 
positions. We use it as a coordinate reference for observations 
without accurate astrometry. 

3. SMA OBSERVATIONS AND DATA REDUCTION 

NGC 1365 was observed using the Submillimeter Array on 
the 4100 m summit of Mauna Kea, Hawaii. The array consists 

7 The other three are NGC 1068, a merger NGC 2146, and NGC 7552 at 
S = -42°. 



of eight 6-meter antennas, cryogenically cooled SIS receivers, 
and digital correla t ors. M ore information on the array is pre- 
sented in Ho et alJ (|2004). Parameters of our observations are 
summarized in Table |2] 

Two nights of observations were made in the 2003-2004 
winter with the array in its compact and extended configu- 
rations. The SIS receivers were tuned to simultaneously de- 
tect three CO lines in the 1 mm band, namely I2 CO(2-l) in 
the upper sideband (USB) and 13 CO(2-l) and C 18 0(2-l) in 
the lower sideband (LSB). The digital correlators have a 2 
GHz bandwidth in each sideband and were configured for a 
spectral resolution of 0.8125 MHz. The southern galaxy was 
observed between the elevations of 15° and 34° under excel- 
lent weather with zenith opacity of 0.05-0.08. The median 
double-sideband system temperature toward the galaxy was 
210 K and 160 K for each night. After flagging bad scans, 
the galaxy data have projected baselines in the range of 6.8 - 
175.0 m and the total integration time of 6.5 hr. 

We made mosaics with a half -power field of view of about 
100" x 50" by observing three positions, namely the galac- 
tic center and two offset positions 25" from the center along 
the position angle of 70°. They were cycled every 12 min- 
utes to give similar u-v coverage for each position. We also 
observed quasars J0132-169 and J0457-234 for gain cali- 
bration, Jupiter, Saturn, Uranus, Callisto, Ganymede, 3C279, 
and J2258-279 for passband calibration, and used Mars and 
Uranus for flux calibration 8 . A small correction, 7% in 
amplitude, is made to our data to compensate for the lack 
of automatic focus tr acking at the time of the observations 
llMatsushita et al.l2006l) . 

The SMA data were reduced in three steps. First, the stan- 
dard passband and gain calibratio ns were made using MIR, 
which is an IDL version of MMA JScoville et alii 9931) . Then 
mosaic images were made from the calibrated visibilities us- 
ing MIRIAD (ISaultetal.lll995l) Finally, the NRAO AIPS 
package (Bridle & Greisen 1994) was used for data analysis. 
In the imaging process, the spectral-line data were binned to 
10 km s" 1 or 20 km s" 1 resolutions and various types of u-v 
weighting were employed. For one set of images, each line 
was given a different u-v weighting (i.e., robust value) in or- 
der to maximize spatial resolution while retaining reasonable 
signal-to-noise ratio. For another, we used exactly the same 
range of u-v radii and the same u-v weighting for different 
lines in order to accurately measure line intensity ratios. Con- 
tinuum data were made by combining channels without line 
emission, and the total bandwidth of the data was 0.7 GHz 
in each sideband. Weak continuum emission was detected at 
about 5a. The continuum was subtracted from the line data in 
the u-v domain. 

We compared our data with previous single-dish observa- 
tions to check consistency. Fig. ^ shows the CO spectra to 
be compared with 25" resolution observations at the galac- 
tic center. They are made from the data cubes used for 
Fig. 12 The asymmetry of the line profile seen in all the 
lines, with stronger blue-shifted emission, is due to an in- 
trinsic asymmetry of gas distribution. The total flux in each 
line is 5.5xl0 3 Jy kms" 1 , 4.1 xlO 2 Jy kms" 1 , and 1.1 xlO 2 
Jy kms" 1 for 12 CO(2-l) , 1 3 CO(2 -l), and C ls O(2-l), re- 
spectively. Sandav ist et al.l ([1995) observed the same posi- 
tion in 12 CO(2-l) with a 25" beam at the Swedish-ESO- 
Submillimeter Telescope (SEST) and obtained a total flux of 

8 Mars was not used in the USB because of the significant 12 CO absorption 
in the sideband. 
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4.9 x 10 3 Jy km s" 1 . Thus the SMA and SEST observations 
agree with each other considering that each probably has an 
uncertainty of ^10%. The SEST line profile has a symmetric 
double peak unlike SMA's, and the small pointing offset im- 
plied by this may be a part of the reasons for the smaller SEST 
flux. We conclude from the comparison that the SMA obser- 
vations detected most of the CO flux in the galactic center. 
We also infer that our 13 CO(2-l) and C ls O(2-l) data sim- 
ilarly detected most of the flux in the observed area because 
they have almost the same u-v coverage as the 12 CO(2-l) data 
and because the u-v weighting used for them are closer to the 
natural weighting. 

In this paper, velocities are with respect to the Local Stan- 
dard of Rest (LSR) and are defined in the radio convention. 
Maps are presented without correcting for the attenuation by 
the mosaicked primary beam, unless otherwise noted. How- 
ever, all the flux and intensity measurements are made after 
the correction. Line spectra are also made from data cubes 
after the correction. 

4. RESULTS 

4.1. Gas Distribution 

The CO maps in Figure|2]show that the emission is mostly 
in the central 4x2 kpc of our 9x5 kpc field of view and that 
brighter emission is on an elliptical ring. There are gas lanes 
on the leading sides of the bar which curve inward near the 
nucleus to form the incomplete circumnuclear ring. This gas 
morphology agrees with that of dust lanes seen in the optical 
(Fig. |3}, except in its western part (i.e., the near side) where 
only the dark lane is conspicuous. The oval ring has a ma- 
jor axis of about 2 kpc (~20"), is made of several distinct gas 
clumps, and is surrounded by weaker emission. Particularly in 
the rare CO isotopes, a pair of condensations stands out near 
the ring's apocenters, where the ring and leading-edge gas 
lanes connect. This feature, resembling the ga s 'twin-p eaks' 
observed in the centers of many barred galaxies ( Kennev et al.l 
119921). is the one previously observed in CO(3-2) and NH3 
( Sandnvist ^T1IT999t lOtt et aUl2005i The CO(2-l) emis- 
sion does not peak at the Seyfert nucleus. The 1.3 mm contin- 
uum is detected at about 5a near the brighter one of the twin 
peaks. 

Radial distribution of the 12 CO integrated intensity, sam- 
pled in concentric circular rings in the galaxy plane and shown 
in Figure^ confirms the molecular ring and the compact size 
of the central gas concentration. The peak at the radius of 
about 5" is due to the ring smeared by the circular sampling. 
The radial profile rapidly declines outwards, falling to 1/10 of 
the peak at a radius of 2 kpc (23")- The extent of this molec- 
ular gas concentration is comparable to that of dust emission 
at 7 and 15 urs\ and that of ionized gas see n at 3 and 6 cm 
iRoussel et alJl200ll iForbes & Norrislll9 98). The warm dust 
and hot gas reflect the active star formation in the region. 

The mass of molecular gas in the galactic center is es- 
timated using a conversion factor from 12 CO integrated 
intensity to H2 column density, Xqo = 0.5 x 10 2l) cm" 2 
(K kms" 1 )" 1 . The adopted coefficient is a middle value 
of the empirical estimates for the center of our Galaxy 
( Sodr oskiet.aU [l995l lOkaet.alJ IT998t iStronget aU l200l . 
and is about a factor of 3-6 smaller than the ones esti- 
mated for the disk of our Galaxy (e.g . ISanders et al.lll984l 
iHunter et al] Il997t iDame et all 1200 lh . The tendency for 
a galactic center, especially a starburst nucleus, to have a 
lower conversion factor than the Ga lactic disk has been sug- 
gested by a number of authors (e.g. lMalonev & Blacklll988t 



iBrvant & Scovilldll996llMeier & Turiierll2001l) . We also as- 
sume that the scaling factor is constant over the region we 
observed, and that the 12 CO(1-0) and 12 CO(2-l) lines have 
the same brightness temperature. Single-dish observations of 
the 2- 1 to 1-0 ratio range from 0.55 to 1.1 for the galactic 
cente r iS andq yis{ et alj 1 98 811 1 995t Papadopoulos & Seaauist 
1998; iCurran et aUl2001l) . The latter approximation makes 
^co(2-i) = ^co(i-O), and hence we omit the transition index. 
For brevity, Xco in units of 10 20 cm" 2 (K km s" 1 )" 1 is written 
asxco m the following; we adopted xco = 0.5. The 12 CO(2-l) 
flux within a 1 kpc radius of the galactic center in the galaxy 
plane is 4.3 x 10 3 Jy km s" 1 . The mass of molecular gas in the 
central 2 kpc is thus estimated to be M mo \(r g < 1 kpc) = 9 x 10 8 
Mq from the above assumptions, the 12 CO(2-l) flux, and a 
factor of 1.36 for the contribution of helium. Molecular mass 
as a function of radius is in Figure^ 

An alternative way to estimate the mass of molecular gas 
uses the CO isotope lines that are probably optically thin. 
The molecular mass in the central 2 kpc is estimated to be 
4 x 10 8 Mq from the C 18 0(2-l) flux of 1.1 x 10 2 Jy km s" 1 
in the region. For this, we used the C ls O abundance of 
[H 2 /C ' 8 Q1 = 1Q 6 - 3 , i.e., rH 7 ./ 12 C01 = 10 4 and r 12 CO/C 18 Q l 
= 200 llHenkel & Mauersberge rl 1 9931 IWilson & Roodi 1 991 . 
and Ar(C 18 0)// c i 80( 2-i) ~ 10 14 6 cm" 2 (K km s" 1 )" 1 suitable 
for warm and moderately dense gas (T^n > 40 K, Nh 2 = 
10 3 -10 5 5 cm" 3 and r (C 18 Q2-l) < 0.1 : Stutzki & Glisten 
H990t lWildetalJ fT992i. This estimate uses the fact that CO 
molecules can be excited only to a limited number of levels 
for the density and, in that circumstance, the level population 
becomes independent of temperature when the gas is warm 
enough. If one instead assumes the LTE, then the mass is 
the same for 20 K gas and larger for warmer gas. The gas 
mass from C 18 is about a factor of 2 smaller than the one 
from 12 CO. This may mean that the adopted Xco is too large. 
However, it is also possible that there is lower density gas 
between or in the outskirts of the molecular clouds, and that 
C 18 can not reliably trace such gas because of difficulty in 
excitation or lower abundance of the molecule in the thin gas. 
For the Galactic center, Dahmenetal. (1998) estimated the 
abundance of this type of thin, invisible, molecular gas to be 
as much as the C 18 0-visible gas. We therefore estimate gas 
mass from 12 CO and the abovementioned Xco through the pa- 
per. For completeness, the 13 CO flux of 3.8 x 10 2 Jy km s" 1 in 
the region gives a M mol of 3 x 10 8 M G for [ 12 CO/ 13 CO] = 40. 
The line is more prone to saturate than C 18 0, and hence the 
mass is more likely an underestimate. The mass-related and 
other key parameters measured or estimated from our data are 
summarized in Table [3] 

The uncertainty in the mass estimates is dominated by that 
in the CO flux-to-mass conversion and may be as large as a 
factor of a few. It is much larger than the ~10% error in the 
CO flux. Though we took the center of our Galaxy as a tem- 
plate for our mass estimate, it is quite possible that the scaling 
relation is offset in the starburst nucleus of NGC 1365. The 
rare CO isotopes give us an independent check, and the es- 
timate of the error. Note that the scaling factor can plausibly 
have local scatter within the galactic center, because gas prop- 
erties are not expected to be uniform throughout the region. It 
is only the lack of observational information sufficient to lo- 
cally determine the gas mass, and the scaling relation, that led 
us to use the constant conversion factor. 

Bearing the abovementioned caveats in mind, we can cal- 
culate the mass and surface density of molecular gas in small 
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regions. Figure |5] shows the distribution of CO integrated in- 
tensity and molecular surface density in the circumnuclear 
region. The maximum integrated intensity of 42 Jy km s" 1 
arcsec" 2 translates to the mass surface density of molecular 
gas S mo i = 9 x 10 2 Mq pc~ 2 in the galaxy plane. The surface 
density is averaged in our beam having 5.1 x 10 4 pc 2 area in 
the plane. The peak is at a compact gas clump 7" north of 
the nucleus. The total mass of molecular gas in the clump is 
^5 x 10 7 M Q . We see later that the gas clump coincides with 
a prominent radio and mid-IR source called D/M4. 

4.2. Kinematics 

The velocity field shown in Fig. ^indicates the presence of 
noncircular motion, most easily seen as the curved (S-shaped) 
velocity contours around the systemic velocity. The noncircu- 
lar motion is, however, rather systematic. The contours show 
a m = 2 symmetry, i.e., they can be rotated by 180° without 
breaking the overall pattern of the velocity field. 

The CO velocity field was fitted with the MPS task GAL to 
derive kinematical parameters; integrated intensity was used 
for data weighting. The large noncircular motion and the 
limited spatial coverage of the CO data preclude us from de- 
termining the inclination and the position angle of the major 
axis. Hence we adopt tho se parameters from H I observations 
of the entire galaxy by Jors ater & van Moorsell (119951) . We 
expect, however, that the systemic velocity and the dynamical 
center of the galaxy are little biased by the symmetric noncir- 
cular motion. The systemic velocity of the galaxy is estimated 
as V sys (radio,LSR) = 1618km s" 1 from the fit, and the dynam- 
ical center of the galaxy is estimated to be a = 03 h 33 m 36?35, 
8 = -36°08'25."8 (J2000). The position coincides with the 
Seyfert nucleus within 0"4, well within the measurement er- 
rors of the dynamical center and the active nucleus. 

The major-axis position-velocity (PV) diagram of 12 CO(2- 
1) emission in Figure [6] shows that the molecular ring is 
around the turnover of the rotation velocity and that no 
anomalous velocity is observed at the active nucleus. The ro- 
tation velocity at the galactocentric radius of 1 kpc (11 "5) is 
approximately 150 km s" 1 , according to the PV diagram. For 
the adopted inclination of 40°, the orbital period at the ra- 
dius and the dynamical mass within the radius are T& yn (R g = 
1 kpc) - 3 x 10 7 yr and M dyn (R g < 1 kpc) - 1 x 10 1() M Q , re- 
spectively. Formal errors in these values can not be properly 
determined without modeling noncircular motion of gas in the 
nonaxisymmetric potential. For a ±20% uncertainty in the 
rotation velocity, which we assume to be reasonable for our 
data, the fractional uncertainties in 7d yn and Md yn are ±20% 
and ±40%, respectively. The gas-to-dynamical mass ratio in 
the central 2 kpc is ~0. 1 . 

4.3. Line Ratios 

The ratio of brightness temperatures is measured between 
12 CO(2-l) and 13 CO(2-l) lines in data cubes of matching 
spatial and velocity resolutions. The line ratio should reflect 
physical parameters of th e molecular gas su ch as density and 
column density (e.g., Mei er & Turnerl20*0 4). and the SMA en- 
ables us to accurately measure it at high spatial resolutions 
through simultaneous observations of the CO lines. Care was 
taken to make the reduction procedures for the two lines as 
close as possible. We used visibilities in the range of u-v radii 
5-128 kA where both lines have data. The same weighting pa- 
rameter, i.e., a robust value of 0.5, and the same velocity bin- 
ning, 20 km s" 1 , were used. The resulting data cubes, which 



had slightly different resolutions, were convolved to the same 
resolution of 3."3 x 2"0. The line brightness temperatures 
were sampled in each velocity channel in a grid of 1 "2 spac- 
ing, which is close to the Nyquist frequency for the spatial 
resolution of the data. The line ratio is calculated only at the 
(a, 8, v) pixels where both lines are detected above 3er. 

Figure shows the scatter diagram of the line intensities 
and the histogram of the 12 CO(2-l) to 13 CO(2-l) intensity 
ratio, Rf~\\ = /( 12 CO 2-l)//( I3 CO 2-1), where /(line) = 

J 7b(line)iiv. The ratio Rfyu * s mostly in the range of 4-14; 
92% of the data are in this range. The mean and the standard 
deviation of Rfyu ^ 8-7 and 2.8, respectively. The former is 

(2— 1 } 

consistent with the single-dish meas urement of Z? 12 /i3 = 10±2 
in the central 21" of the galaxy ( Papadon oulos & Se aauist 
1998). The latter is an upper limit of the intrinsic variation 
of the ratio because of the noise in the data. Small number of 
pixels having extreme ratios are likely due to noise. In par- 

(2— 1 ) 

ticular, pixels having /?j 2 / 13 < 2, which are less than 2% of 
the data, are almost certainly noise judging from their mostly 
isolated locations in the channel maps. 

(2—1) 

The spatial distribution of the ratio ^ 12 /i3 * s shown in Fig- 
ure [8] For this map, line ratios are averaged at each position 
if there are multiple velocity channels with a valid ratio. The 
overall trend in the map is that the ratio is relatively high at the 
ridges of molecular gas. Specifically, higher ratios are seen on 
the southern molecular ridge at the southern dust lane in the 
bar, on the southeastern half of the molecular ring, and on the 
northern molecular ridge in the bar. Lower ratios are mainly 
observed in the northwestern half of the circumnuclear ring 
and inside the ring. The 12 CO(2-l) to C 18 0(2-l) ratio has the 
same trend, as seen from the comparison of Figures|2](a) and 
|2](e). Another way to describe the spatial distribution of the 
line ratio is that there is a trend of the ratio being higher on 
the far (SE) side than the near (NW) side in the circumnuclear 
disk of 1 kpc radius (see Figs. [7Jand[8j. The mean and stan- 

(2—1) 

dard deviation of ^ 12/ / I3 are (mean ± a) = 7.9 ± 2.2 for the 
near side and 10.0 ±2.6 for the far side. A similar trend of 
higher ^ 12 /i3 in the far side has been noticed in the circumnu- 
clear regions of IC 342 and NGC 253 (iMeier & Turnerl2001l 
Sakamoto et al. 2006a), but the illumination model proposed 
for them may not apply to NGC 1365. This is because the 
starburst ring coincident with the molecular ring can not pref- 
erentially illuminate the side of clouds facing us, nor is there 
a radial gradient of ^12/13 expected from AGN illumination. 
We further discuss the ratio in ^6] 

4.4. Comparison of CO with Radio and mid-IR sources 

We show in Figure[9]the locations and conventional names 
of the compact radio sources, luminous mid-IR sources, and 
three of the optically identified super clusters in the circumnu- 
clear region. The radio sources have alphabet names, and the 
mid-IR sources are called M#, where # is a number. As seen 
in the figure, some of the radio and mid-IR sources coincide 
with each other, and are thought to be the sam e objects, likely 
young massive clust ers embedded in dust (Sandavi st et alJ 
n995TlCTalliano et al.l l2~005). In contrast, the mid-IR sources 
M2 and M3 do not have a radio counterpart but coincide with 
optically iden tified clusters called SSC 3 and SSC 6, respec- 
tively, which Rristen et al. (1997) found to be the brightest 
optical sources in the region except the nucleus. Their mid-IR 
colors suggest less dust emission than toward other sources 
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JGalliano et al.l2005T) . consistent with their optical detection. 

Our CO data are compared with the radio and the mid-IR 
sources in Figure [H)] We have produced for the comparison 
peak brightness temperature maps in addition to the CO in- 
tegrated intensity maps. Since the peaks in these maps often 
coincide, we collectively call the peaks CO hotspots. In the 
maps, radio sources are plotted as squares at average positions 
of previous observations 9 , while mid-IR sources are marked 
as crosses. 

The comparison shows that some of the radio peaks and 
mid-IR sources coincide with CO hotspots, i.e., peaks of CO 
integrated intensity or CO brightness temperature or both. 
Specifically, the sources D/M4, E/M5, G/M6, and H have a 
local maximum of 12 CO integrated intensity within 1". A 
local peak in 13 CO integrated intensity is also seen within 
1" of D, E, and G. Each of the four sources also coincides 
with a local maximum of peak brightness temperature in I2 CO 
and 13 CO, and the source D is at the maximum of the C ls O 
brightness temperature map. Though not coincident with an 
isolated peak, the radio source F is within 1" of a ridge of 
12 CO integrated intensity, and is on a sharp ridge of 12 CO peak 
brightness temperature. In contrast, the sources M2/SSC6 and 
M3/SSC3 that are thought to be less embedded in dust are in- 
deed in the areas of less molecular gas. The radio sources 
A and J do not have a nearby peak in the CO maps. Table 
0] summarizes the associations of peaks in CO maps to the 
sources in other wavelengths. In general, the correlation to 
the sources appears better in peak CO temperature than in the 
line integrated intensity. The weak 1.3 mm continuum does 
not coincide with any of the radio and mid-IR sources. 

The CO spectra measured at the hotspots are predominantly 
single-peaked, and the peaks of different isotopic lines are at 
the same velocity for each hotspot, as shown in Figure ^2 
This suggests that each hotspot is a single subregion in each 
gas complex. The properties of CO emission measured at the 
CO hotspots are summarized in Table [5] 

4.5. Seyfert Nucleus 

The active nucleus is in the CO deficient area inside the 
molecular ring, and no CO peak is at the nucleus. The inte- 
grated intensity and peak brightness temperature of 12 CO(2- 
1) at the AGN position are 9.3 Jy km s" 1 arcsec" 2 and 2.0 
K, respectively, in our 2. "6 x 1"8 resolution data. The sur- 
face density of molecular gas in the galaxy plane is 2 x 10 2 
M Q pc~ 2 in the central 200 pc for the adopted xqq. The 
mean total extinction in our line of sight is ^12 mag (i.e., 
A/h, = 1 .2 x 10 22 cm -2 ). The column density and extinction in 
front of the AGN would be half of these if the nucleus were at 
the mid-plane of a uniform gas layer. The mass of molecular 
gas within 100 pc of the nucleus is 6 x 10 6 M . The supposed 
molecular torus confining the conical plasma outflow should 
be in t his region, considering th e opening angle of the flow 
( ~50°: iHielm & Lindbladll 19961) and the typical scale height 
of the molecular gas layer (< 100 pc) in the centers of disk 
galaxies. Note that these mass-related quantities inherit the 
uncertainty of the conversion factor, and that the factor quite 
possibly changes between the starburst ring and the vicinity 
of the AGN. 

9 T he r adio source H was g iven positions 1"4 apart by Sand qvist et al] 
fH%)5l) and lStevens eTal] tlTOI) . We use the average position and associate 
the source with M8, which is O'.'l from the radio position. There remains 
a possibility, however, that there are two radio sources and M7 and M8 are 
associated with each of them. 



We detect no kinematical anomaly at the nucleus, such as 
high velocity gas suggestive of fast rotation around the AGN 
or a kinematical sign of a molecular outflow that may be en- 
trained by the plasma outflow. Our spatial resolution, how- 
ever, is probably insufficient to see gas motion dominated by 
the central object. The full-width at zero intensity of 12 CO 
line is 270 km s -1 at the nucleus. This implies a dynamical 
mass of Mdy n (Rg < 0.1 kpc) ~ 1 x 10 9 M Q for gas in circular 
motion in the galaxy plane. This is an order of magnitude 
larger than the mass of the central engine i nferred from the 
K magnitude of the galaxy by Risaliti et al. (2005). Judging 
from the PV diagram (Fig. |6), a factor of >3 higher resolu- 
tion is likely needed before one starts to see the gas motion 
influenced by the expected black hole. 

The properties of molecular gas around the AGN can be 
constrained, though weakly, from the low brightness temper- 
ature of 12 CO. The physical temperature of molecular gas in 
the central 200 pc would be 6 K if the gas is uniformly dis- 
tributed in the region and if I2 CO is optically thick and in 
LTE. Judging from the low temperature, it is more likely that 
one or more of these conditions is not satisfied; i.e., the gas 
has clumps, a 12 CO(2-l) opacity < 1, subthermal excitation 
of CO, or a combination of these. The clumpy distribution of 
the ISM could reduce the extinction toward the AGN, to be 
more consistent with the optically visible nucleus. 

5. BAR-DRIVEN GAS DYNAMICS 

The central concentration of molecular gas, with an ex- 
tent of a few kpc and a mass of ^10 9 M©, is almost an ex- 
pected feature for NGC 1365 considering the long and promi- 
nent bar in the galaxy. A stellar bar is known to drain the 
gas in the region it swee ps toward the g alactic center (e.g., 
Matsuda & Nelsonl U977t ISakamoto et a"l1ll999bt LSheth et al] 
2005). The bar in NGC 1365 has apparently worked as a pow- 
erful gas funnel. 

The overall distribution of molecular gas in the central con- 
centration can be understood in the framework of the bar- 
driven gas dynamics, in which the dominant pattern of gas 
flow changes between the so-called jti-type and X2-type in the 
circumnuclear region. The x\ and X2 are families of stable, 
prograde, and periodic orbits for st ars in the nonaxisymmet - 
ric potential of a galaxy bar (see [Binnev & Tremaine 1987). 
The former orbits are elongated along the bar and fill most 
of it, and the latter, which can exist closer to the galactic 
center under certain conditions, are elongated perpendicular 
to the bar. Gas follows similar orbits or streamlines, but 
because gas is collisional unlike stars the angle between x\ 
and X2 for gas can be oblique and less than 9 0° measured in 
the direction of galaxy rotation (Wada fl994l). As has been 
mo deled specifically for t his galaxy by iTeuben et alJ (119861) 
and Lindblad et al. ( 1996) and for generic barred galaxies by 
many authors, if there exist xi orbits, gas can form a circum- 
nuclear ring or spiral pattern connected to a pair of straighter 
gas lanes in the bar. What we see in NGC 1365 matches the 
model gas morphology and what has been seen in many other 
barred galaxies. The leading-edge gas lanes on the bar are part 
of the jci -like orbits and the gas ring, whose major axis makes 
an oblique angle less than 90° with the bar, corresponds to the 
X2-like orbits. 

Kinematic ally, there are two pieces of evidence that sup- 
port the gas dynamical model mentioned above. First, the 
S sha ped isovelocity contours around the systemic velocity 
( 34.2> are explicable by the transition of dominant orbits be- 
tween jci -like outside the ring and X2-like around the ring as 
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seen in the models with the x\—X2 transi tion, or an inner Lind- 
blad resonance (ILR), in Fig. 4 o f lTeubenlTaTl fiWh and 
Fig. 20b of Lindblad et all ( I1996T) . The noncircular m otion 
can b e caused purely by the transition jLindblad & Lin dblad 
1994), but gas shocks may be also contributing. Secondly, 
as shown in Fig. El the double-peaked line profiles near 
an apocenter of the circumnuclear ring are consistent with 
oval gas motion along the ring, whose major axis makes an 
oblique angle to the bar. This again suggests that X2-like 
gas stream exists on or slightly inside the molecular ring. 
We note that these kinematical and morphological clues for 
the presence of a x\-x 2 transition are probably more robust 
than an analysis of the rotation curve by using the resonance 
condition of 2(0 - fib) = where Q, f^, and k are angu- 
lar frequency, the pattern speed of the bar, and epicyclic fre- 
quency, respectively. The resonance condition is known to be 
a poor indicator for the presence of x? o rbits in a strong bar 
dContopoulos & Papavannopoulos 1983). 

6. GAS PROPERTIES 
6.1. Average Properties in the Galactic Center 

The observed 12 CO-to- I3 CO intensity ratio, Rf' 1 ^ = 8.7 ± 
2.8, is normal for the centers of local spiral galaxies including 
the ones with moderate starburst (log(Lg_iooo ^hi/Lq) < 10.5). 
For example, the ratio is about 10 at the centers of our G alaxy, 
IC342, NGC 253, and M82 (ISawada et alJl2001tlMeier et all 
l2000HMao et al.l2000USakamoto et al.l2006al) . The ratio sep- 
arates NGC 1365 from the more infrared luminous mergers 
(log(L8_iooo (uh/Lq) > 1 1 .5), which are known to show higher 
ratios of Rf;^\ > 20 l|C^.oli et al.l[l992t lAalto et al.l lT995l 

Glenn & Hunt erl200lHSakamoto et all2006bl) . 

There are three basic circumstances under which the I2 CO- 
to- 13 CO intensity ratio is about 10, although some of their 
intermediate situations also produce the same ratio. The three 
correspond to the three parts of the formal expression of the 
ratio: 

_ (l-e-^) [J(T exA2 )-J(T CMB )] f bA2 
12/13 (l-e-^) [/(^WCTcmb)] / b>13 ' 

where r, T ex , and /b are optical depth, excitation tempera- 
ture, and beam filling factor, respectively, and J{T) is the ef- 
fective radiation temperature (hv /k)/[exp(hv /kT)- 1] for the 
line frequency denoted as v. First, if the CO is thermalized 
at least up to the J=2 level and if the two isotopic lines have 
almost the same beam-filling factors, then the ratio is deter- 
mined by the optical depths. For example, the ratio of 8.7 
corresponds to ri3 CO (2-i) = 0.12 and 13 CO column density of 
5 x 10 15 cm -2 (km s" 1 for a 50 K gas in LTE. The molecular 
gas needs to be turbulent, with a characteristic velocity gradi- 
ent of dv/dr ~ 15 km s" 1 pc" 1 , to have a hydrogen density ex- 
ceeding n H , ~ 10 4 cm" 3 to maintain the excitation. Secondly, 
the ratio of ~10 can also be observed when the 13 CO excita- 
tion is highly subthermal. In t he non-LTE analysis using the 
large velocity gradient model dGoldreich & Kwanlll974l) . the 
mean ratio of 8.7 requires a gas density of n^ ~ 3 x 10 2 cm" 3 
for the abundances of [H 2 / 12 CO] = 10 4 and [ 12 CO/ 13 CO] = 40 
and a velocity gradient in clouds of 1 km s" 1 pc" 1 . Gas tem- 
peratures in the range of 10 - 300 K produce the same result. 
In these conditions, the excitation to J=2 is subthermal with 
13 CO having lower excitation temperature than 12 CO, and the 
13 CO(2-l) opacity is 0.5-0.9. Finally, different beam filling 
factors of the two lines can affect the ratio if molecular clouds 



in the galactic center have a core-envelope structure as has 
been observe d in the Galactic center and sug gested in lumi- 
nous mergers (Huttemeis ter et all 19951 lAalto et all!995l) . In 
this case, 12 CO traces lower density envelopes while 13 CO is 
emitted from higher density cores. 

Our mass estimate of molecular gas assumes the nearly 
thermalized excitation of the first case rather than the very 
subthermal one in the second. In the latter, C I8 0(2-1) exci- 
tation temperature is about a tenth of the thermal gas temper- 
ature for a 50 K gas. Even though we already allowed for 
the contribution of thin gas, our mass estimates would need 
to be doubled if most of the C 18 and 13 CO emission is from 
very poorly excited molecules. Large turbulence in molecular 
gas, which separates the thermalized case from the subther- 
mal one, has been suggested t o be a cha racteristic feature of 
the Galactic center (Sawadaetal. 2001). The turbulence is 
measured with a velocity gradient in the escape probability 
approximation and is larger than 10 km s" 1 pc" 1 . We expect 
the starburst nucleus of NGC 1365 has no less turbulence than 
the Galactic center, and hence keep our mass estimates. The 
very low excitation of C ls O and CO can not be ruled out, 
however, without further observations. 

6.2. Spatial Variation 

(2— 1 } 

The variation of n 12 /i3 within the circumnuclear region, by 
about a factor of 3 at our resolution, allows multiple interpre- 
tations corresponding to the three cases mentioned above. In 
the thermalized case, a factor of 3 lower ratio can mean a gas 
column density per velocity (A^h, / dv) increased by the same 
factor, or a gas temperature reduced by about a factor of 2. 
In the subthermal case, a higher gas density reduces the ratio. 
In the filling-factor model, molecular clouds in lower R\2/n 
regions have smaller envelopes with respect to their cores. 

(2—1) 

Our data show the lack of an organized pattern in the R\ 2 / l} 
ratio, except for the slight difference between the near and far 
sides. This makes it difficult to assess the gas properties by 
combining the line ratio and other information. For example, 
it is possible that molecular gas is relatively dens e in the twi n 
peaks owing to the convergence of gas flow (lOtt et all 2005). 
One of them, the gas complex containing radio sources D, 
E, and G, shows a lower ratio of ~8 than other bright 12 CO 
complexes, which have ratios > 10 (Fig. [8}. This is consis- 
tent with the higher density, which also makes N^Jdv higher 
for a constant turbulence (dv/dr). However, there are regions 
having even lower ratios without obvious dynamical reasons. 
Quite possibly, multiple parameters change at the same time 
for multiple reasons, making any large-scale ratio pattern less 
clear. 

On smaller scales, our data also show no distinctive corre- 

(2—1) 

lation between R l2 / l3 and the radio, optical, and IR sources. 
This likely means that the effect of these sources on the sur- 
rounding molecular gas is either small in terms of the ratio or 
is heavily diluted in our 200 pc beam. The former is possible, 
for example, if higher gas temperature around a cluster, which 
could reduce the optical depths of the CO lines for a constant 
column density, is offset by the gas concentration (i.e., larger 
column density) around the clusters. The ratio Ru/13, which 
traces the 13 CO optical depth in the thermalized case, would 
then be comparable toward the clusters and the surrounding 
area. Finally, it is noteworthy that the high ratio of > 20 ob- 
served in luminous mergers in kpc beams is not observed in 
any part of the starburst region. 
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7. CIRCUMNUCLEAR STAR FORMATION 
7.1. Star Formation Law 

The circumnuclear starbur st in NGC 1365 follows the two 
empirical scaling relations of lKennicuttl (I 1 9981 hereafter K98) 
for starburst galaxies. Taking 1 kpc as the radius of the star- 
burst region, and attributing half of the far-IR luminosity of 
the galaxy to the starburst (see the star formation rate in 
the region is 9 Mq yr" 1 and its surface density is Esfr = 3 Mq 
yr" 1 kpc" 2 for the star formation rate -to- far-IR luminosity re- 
lation of K98. The mean surface density of H2 in the region is 
Eh 2 =1 x 10 3 (x CO /2.8) M Q pc" 2 ; K98 adopted x co = 2.8. Fig. 
H3lshows that NGC 1365 with these surface densities follows 
the correlation, or the Schmidt law, seen among circumnu- 
clear starbursts. The galaxy also follows the same correla- 
tion for a combined sample of starbursts and normal galaxy 
disks (see Fig. 6 in K98). In addition, NGC 1365 follows 
another correlation between gas surface density normalized 
by the dynamical time and the star formation surface density. 
The former is S H2 /r dyn =4 x 10 1 (x C o/2.8) Mq yr -1 kpc" 2 for 
the starburst region having the 7d yn calculated in 34.21 This 
and the abovementioned Esfr place NGC 1365 on the tight 
correlation between the two quantities seen in Fig. 7 of K98, 
with the abscissa and ordinate being 3.6 and 0.5, respectively, 
in the figure's units. 

The fact that the circumnuclear starburst in NGC 1365 fol- 
lows these empirical correlations is not affected by the choice 
of the conversion relations used to derive H2 mass and star 
formation rate, because we use the same conversion factors 
as in K98. Removing the conversion relations and stated in 
observed quantities, the circumnuclear starburst in the cen- 
tral 2 kpc of NGC 1365 follows the correlation, seen among 
starbursts and galactic disks, between far-IR luminosity av- 
eraged over the star forming region and CO ( J= 1 — equiva- 
lent) luminosity averaged in the same area. The second cor- 
relation that NGC 1365 follows can be also stated only us- 
ing observed quantities. The factor of ^6 difference in the 
CO to H2 conversion factor between this paper and K98 only 
affects such quantities as gas consumption time and gas-to- 
dynamical mass ratio. The time in which the current amount 
of gas would be consumed by the star formation at the current 
rate is 1 x 10 8 yr in the 2 kpc-diameter starburst region for 
the conversion factor we adopted, and it is 6 times long er for 
the disk X C o used in K98. Note that iKennicuti (119981) used 
the disk conversion factor for his sample only for the sake 
of simplicity, being fully aware that the factor can change in 
starburst environment. 

7.2. Reasons for the Large Luminosity 

The reasons for the starburst in NGC 1365 being one of 
the most luminous in non-merging galaxies while following 
the normal star formation laws are the high gas surface den- 
sity in the starburst and the large size of the starburst region. 
The mean gas surface density in the starburst is among the 
highest observed in non-merging starburst galaxies, as seen 
in Fig. 1131 where the galaxy has the fourth highest Sh,, 
only after NGC 3079, M82, and NGC 253, among 17 non- 
merging starbursts. The central gas-surface density of NGC 
1365 also stands out when compa red to those in t he NRO- 
OVRO survey of nearby galaxies (Sakamot o et alJlT999albl) . 
The high-resolution survey gives gas surface densities in the 
central kiloparsec of 20 spirals, including some starbursts but 
no merger. The gas surface density in the central kpc of NGC 
1365 is Z mol (R g < 500 pc) = 4 x 10 2 (x co /0.5) Mq pc" 2 and 



is twice larger than the largest in the survey if the same xco 
is used for all galaxies. Meanwhile, the extent of the circum- 
nuclear starburst in NGC 1365, ^2 kpc, i s among the largest 
compared t o those in K98 and s urveys bv lJoeee et alJ (2005. 
Fig. 5) and lKnapen et alJ d2006l Fig. 1). Obviously, the high 
gas surface density in the large are a nee ds a large amount of 
gas; M mo i(r g < 1 kpc) ~ 10 9 Mq ( 34. H . For this accumula- 
tion, the large size of the galaxy (L B » 6 x 10 10 L Q « 3L£, 
D25 = 59 kpc; de Vaucouleurs et all 1 19911) and the long bar 
('bar ~ 18 kpc) were probably crucial through the bar-driven 
transport of disk gas. Thus the combination could be regarded 
as a deeper cause for the luminous starburst. 

The luminous starburst in NGC 1365 contains subregions 
comparable to the entire starbursts in M82 and NGC 253. 
Most notably, the northeastern edge of the circumnuclear ring 
has Sh 2 ~3x 10 3 (xco/2.8)M q pc" 2 , or the same gas surface 
density as in those two starbursts for the same extent of ~400 
pc and for the same xco- The radio and mid-IR sources in the 
region suggest enhanced Ssfr there. With a few such subre- 
gions, it is reasonable that NGC 1365's starburst is a factor 
of 2-4 more luminous than those in M82 and NGC 253. Al- 
though the empirical star formation laws in the previous sec- 
tion and the scaling observation like this leave out the physics 
of star formation on smaller scales, they are still useful to ob- 
tain an overview — star formation in this luminous starburst 
is not discrete in nature from those in less luminous starbursts. 

It remains to be found out why the luminous starburst 
is seen now even though it can last only for ^10 8 yr. 
Rov & Walsh ( 1997) suggested that the bar was young (age 
~1 Gyr) and started the gas transport in the recent past on 
the basis of a bend in the radial gradient of O/H abundance, 
but the existence of the bend has been disputed by Pilvugin 
( 120031) . The galaxy is a member of the Fornax cluster, but 
lacks an obvious interacting partn er for strong gravitational 
disturbance. H I imaging by Jtirsater & van Moorsel (IT9951) 
does not reveal a sign of recent galaxy interaction either, 
though the outer H I disk shows a hint of interaction with the 
intergalactic medium. It may be that the starburst has been 
turning on and off without external triggers, per haps owing 
to a t hreshold in the star formation mechanism ( Jogee et al. 
2005). If so, the star formation law must switch between the 
one for the starburst state (e.g., Fig. 1 1 31 and the other for 
non-starburst. 

7.3. CO Hotspots and Super Clusters 

We saw in §!4.4l that the molecular ring contains several CO 
hotspots and that three of them are spatially associated with 
the radio/mid-IR sources D/M4, E/M5, and G/M6 that are 
suggested to be sup er clusters of ~10 6 M q, ^10 9 Lq and 
ages of 3-6 Myr bv lGalliano etaD J2005I) . Another source, 
H/M8, also coincides with a CO hotspot, but it is fainter and 
does not have age and mass estimates. It is possible that each 
of the radio/mid-IR sources is a group of clusters rather than 
a single cl uster, since their sizes are only constrained t o be < 
20-60 pc (SandQvi stetai]fl995l iGaihano et alJ l2005) while 
super clusters have sizes of < 10 pc. We refer each source as 
a super cluster for simplicity, but 'a cluster' in the following 
can be 'a cluster complex'. 

The empirical star formation laws that we found applica- 
ble to the circumnuclear region as a whole allow the amount 
of star formation needed for the clusters. A star forma- 
tion rate of Ssfr ~ 20 Mq yr" 1 kpc" 2 is expected from the 
star formation laws in K98 for the typical CO integrated in- 
tensity on the circumnuclear ring, 7co2-i ~ 150 Jy km s" 1 
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beam equivalent to an inclination corrected surface density 
of Sh 2 ~ 3 x 10 3 (x C o/2.8) M Q pc -2 . Thus, if the star forma- 
tion laws apply to 100 pc scales, 10 6 Mq of stars will form 
every 5 Myr in every 100 pc by 100 pc area in the molecular 
ring. The empirical laws, however, do not tell whether and 
where a compact and massive cluster forms. The spatial asso- 
ciation of the putative ^10 6 Mq clusters and the CO hotspots 
may be a valuable observational clue to address these points. 

The significance of the hotspot-cluster association is that 
the CO hotspots were very likely the sites of cluster formation 
in the recent past for two reasons. First, the clusters' relative 
motion with respect to the ambient gas in their ~5 Myr lives 
would be < 50 pc if we assume the clusters' relative velocities 
against the ambient gas to be ~ 10 km s , which is a typical 
velocity dispersion of molecular clouds in galaxies. Accord- 
ing to this estimate, the CO hotspots seen in our ~200 pc res- 
olution images are composed of gas that was in the immediate 
vicinity of the formation sites of the super clusters. Second, 
given the non-trivial configuration of the three hotspot-cluster 
pairs, it is unlikely that the three clusters formed elsewhere 
and have migrated to the CO hotspots, unless the clusters 
themselves affected their ambient gas to create the hotspots. 
This argument reinforces the first one, in which the drift ve- 
locity might be larger than assumed because the three clusters 
are in the region with two types of orbits and because only 
gas is subject to hydrodynamical forces. Regarding the possi- 
bility that the clusters caused the hotspots rather than the gas 
hotspots formed the clusters, we favor the latter interpretation 
though our data do not rule out the former. Our preference 
is because three hotspots are visible in 13 CO and one even 
in C ls O. Elevation of gas temperature, the most likely effect 
of a young cluster on its ambient gas, does not increase the 
intensity of an optically thin CO line per unit mass of molec- 
ular gas if the isot ope CO is well excited, which we favorably 
discussed in 36. II 

In our more favored model, the CO hotspots are peaks of 
gas surface density. The surface density is either that for the 
entire velocity or for a 10 or 20 km s -1 range depending on 
whether the hotspot is in the integrated map or in the peak 
temperature map. Gas temperature can also be higher in the 
hotspots seen in I2 CO because the brightness temperature of 
the line is the gas temperature multiplied by the beam filling 
factor. The gas surface densities in the hotspots are S mo i ~ 10 3 
Mq pc~ 2 when averaged over 200 pc scale (see Fig. [3}. The 
CO hotspot-cluster associations suggest that the super clus- 
ters formed in such a high surface-density environment. If 
the high surface density of molecular gas is one of the re- 
quirements for the SSC formation, then SSC formation pre- 
dominantly in starbursts can be a consequence of the high gas 
surface densities in starburst regions. 

In the alternative interpretation, the CO hotspots are largely 
due to higher gas temperature there, plausibly owing to the lu- 
minous super clusters. For optically thick I2 CO, this straight- 
forwardly explains the hotspots. For 13 CO and C I8 0, this in- 
terpretation is possible only when the gas excitation is sub- 
thermal and the isotope lines h ave moderate (~1) optical 
depths as in the second case in 36.11 In such a situation, a 
higher gas temperature can increase the excitation temper- 
ature of the subthermally-excited CO without significantly 
reducing the line optical depth, and hence can increase the 
line brightness temperature. The gas surface densities at the 
hotspo ts in this case would be higher than in the thermalized 
case ( 36. 1> . even though the hotspots may not be local peaks 



of gas surface density. 

7.3.1. Similar Sources in Other Galaxies 

There have been a small number of cases, in other galaxies, 
that are similar to the CO hotspot-cluster associations. Table 
[6]compares them with those in NGC 1365. These cases show 
that at least some of the most massive and youngest clusters 
in starbursts are accompanied by CO hotspots detectable at 
^100 pc resolutions. 

One case is in the Antennae merger (NGC 4038/4039). The 
brightest 15 /im source in the system, located in the overlap 
region of the colliding galax ies, is associated with a molecular 
comp lex called SGMC4-5 JWilson et alJl2000t iMirabel et alJ 
1998). Their parameters are within a factor of a few from 
those of the CO hotspot-cluster associations in NGC 1365 
considering the larger CO beam and hence more beam di- 
lution in the Antennae. The source has been suggested to 
be a young (~4 Myr), massive (~10 7 Mq), and compact 
(^haif light ~ 30 pc) cluster or cluster complex with thousands 
of O stars and a total luminosity of ~10 9 Lq iGilbert et all 
120051 iNeff& Ulves'tadl l200Cl iHaas et alJ 12000) . Since the 
cluster is in the molecular complex but not at its CO peak, 
the CO hots pot is most l i kely a peak of gas surface density. 
Statistically, Zhang et al. (2001) showed that young clusters 
(of ages ^5 Myr, including the one mentioned above) in this 
galaxy are better associated with mid/far-IR, CO, and radio 
sources than older clusters. 

Another case is in the barred starburst spiral M83. The 
brightest CO peak in the galaxy's circumnuclear region, in 
both integrated intensity and brightn ess temperature, coin- 
cides with a radio and mid-IR p eak (Sakamoto et al] 120041: 
Turne r & Holll994t lTelesco et al. 1993). The radio and mid- 
IR emission is due to star formation, though there may 
be contribution fro m a radio sup ernova to the radio peak 
burner & Holll994UTelesco et alJl!993l) . The CO hotspot is 
on a circumnuclear gas ring and is where xi-like and X2-like 
gas streamlines converge. The hotspot location with respect 
to the bar and the circumnuclear gas ring is therefore identical 
to that in NGC 1365. 

7.4. Cluster and Hotspot Evolution 

It has been proposed that gravitational collapse in an 
ILR gas ring, where a bar accumulates disk gas, forms 
gas clump s and subsequently causes intense sta r formation 
in them jwada & Habelll992t lElmegreed ll994). The CO 
hotspots with embedded young super clusters may be such 
gas clumps. If the molecular hotspots are high surface-density 
regions of molecular gas, possibly caused by gravitational in- 
stability, then the formation of super clusters from such giant 
gas concentrations of ^100 pc and >10 7 Mq scales is in ac- 
cordance with the supergiant-cloud model for the formation 
of super clusters (IWilson et alJl2003l and references therein). 
The CO hotspot-cluster association, however, is also compati- 
ble with another model for SSC formation where a super clus- 
ter forms from a moderate-size mole cular cloud (5, 1 Q 6 M ( t i ) by 
converting most of the gas into stars l lKeto et al. 2005, and ref- 
erences therein). Such clouds quite possibly exist in the high 
surface-density environment of the CO hotspots and could 
have had such triggers as cloud collision and compression 
by nearby H II regions or supernova remnants. High surface- 
density clouds with a few 10 pc sizes and velocity widths of 
the order of 10 km s -1 could cause the CO hotspots in the peak 
temperature maps. 
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The SSCs in the hotspots are expected to disperse ambi- 
ent gas as they age while blowing stellar winds and hav- 
ing supernovae. A super cluster may cause an expanding 
bubbl e or even an ou tflow of the ISM from the disk (e.g., 
Saka moto et al.ll2006albl in our SMA survey). Such features 
may quite possibly be hidden within the CO hotspots. The 
optically-visible super clusters that are not associated with 
CO hotspots, most notably M2/SSC6 and M3/SSC3, may be 
older super clusters that have already dispelled their dust co- 
coons in which they formed. Alternatively, they could have 
formed from isolated molecular clouds and have mostly con- 
sumed them at the time of formation. 

8. SUMMARY 

We have observed 1.3 mm CO lines and continuum in the 
center of the barred spiral galaxy NGC 1365 at high spa- 
tial resolutions. The circumnuclear starburst in the galaxy is 
among the most luminous and intense in non-merging galax- 
ies. Our main findings regarding the properties of molecular 
gas and their relation to the starburst and to the Seyfert nu- 
cleus are as follows. 

1 . There is an oval ring of molecular gas around the nu- 
cleus, having a ~2 kpc extent and twin-peak enhance- 
ments at its apocenters, and connected to gas lanes on 
the leading side of the stellar bar. The gas morphol- 
ogy and kinematics in the region are consistent with the 
model that the ring is at the transition region between 
X[- and X2-like gas stremlines in the bar. The structure 
is typical of barred galaxies. 

2. The molecular ring coincides with the starburst ring 
seen in radio emission, suggesting that the gas distri- 
bution largely determines the extent and the shape of 
the starburst. Moreover, the intensity of the starburst 
averaged over the circumnuclear region is what is ex- 
pected from the molecular gas surf ace density of th e 
region and the star formation laws in Kennicutt (1998). 
A plausible synopsis of the galaxy's starburst is as fol- 
lows: The long (~48 kpc) stellar bar funneled a large 
amount of molecular gas (~10 9 M Q ) toward the galac- 
tic center, and the resulting high gas surface density 
(~10 3 Mq pc" 2 ) in the relatively large circumnuclear 
ring (r~l kpc) resulted in the starburst (~3 M Q yr" 1 
kpc" 2 ) of large total luminosity (~10 n L Q ) following 
the star formation laws. — A remaining question is 
why the transient starburst (gas consumption timescale 
^10 8 yr) is seen now. 

3. There are compact (< 200 pc) peaks of CO integrated 
intensity and brightness temperature, which we call 
CO hotspots, on the circumnuclear ring. Some of the 
hotspots coincide with the radio and mid-IR sources 
that have been identified as dust-enshrouded super star 
clusters. The association of the CO hotspots and super 
clusters allows two interpretations. The first one is that 
the most luminous and presumably massive super clus- 
ters were formed in regions with the highest gas surface 
densities (~10 3 M Q pc~ 2 in 200 pc scale). The second 
one is that the luminous embedded clusters heated their 
ambient gas to form the hotspot. The former, which 
is more conventional in terms of the interpretation of 
CO intensities at > 100 pc scales, implies a relation be- 
tween the cluster formation and gas complexes already 
identifiable at a 200 pc resolution. 



4. The active nucleus coincides with the dynamical center 
of the galaxy within 30 pc, according to the CO ve- 
locity field. Neither integrated intensity nor brightness 
temperature of CO peaks at the nucleus. We found no 
anomaly in molecular gas that may be related to the 
AGN. 



We thank the SMA staff for their help prior to and dur- 
ing our observations, Jun-Ichi Morino and Andrew Baker for 
stimulating discussions, Raffaella Morganti for kindly provid- 
ing her radio data, and the anonymous referee for his/her help- 
ful comments. This research made use of NED, ADS, and the 
HST archive, in addition to the SMA; we are grateful to their 
operating, funding, and contributing people and agencies. 
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TABLE 1 
NGC 1365 PARAMETERS 



parameter 


value 


ref. 


distance [Mpc] 


17.95 


1 


scale 


l"=86pc, 11." 5=1 kpc 




inclination [deg.] 


40 


2 


major axis P.A. [deg.] 


220 


2 


bar P.A. [deg] 


90 


3,4,5 


nucleus" R.A. (J2000) 


03 h 33 m 36?38 


5,6,7 


Dec. (J2000) 


-36°08'25y7 


5,6,7 


l°gU-8-1000 fim/LQ) 


11.00 


8 


Hubble type 


SB(s)b 


9 


References. — 


1. JFreedman et all 12001ft: 2. 



3. 



4. IReg an & Elmegreen 1 19971; 5. Measured in th e 
data of IJarrett et all 120031) : 6 ISandavist et all fTpS^ ; 
7. IStevensetall <1999ftT 8. ISanders et all 120031) : 9. 
Ide VaucouTeursetall 11991ft 

a The adopted mean position is within 0f'45 of the posi- 
tions in three wavelengths. 



TABLE 2 
SMA OBSERVATION PARAMETERS 



parameter 

date (UT) 

array configuration 

No. of antennas 

T225 a 

<7/ sys (DSB)) b [K] 
S„(J0132-169) C [Jy] 
S„(J0457-234) C [Jy] 
baseline length [m] 
integration time d [hr] 
center frequency, L/U [GHz] 
bandwidth [GHz] 
spectral resolution [MHz] 
mosaic center (J2000) 

primary beam c [arcsec] 
mosaic field of view [arcsec 2 ] 
rms noise f [mJy beam -1 ] 



value 



2003-12-23 2004-02-02 



compact 


extended 


6 


7 


0.08 


0.05 


210 


160 


1.38/1.35 


1.28/1.20 


1.28/1.20 


0.77/0.75 


6.8-64.9 


11.4-175.0 


3.8 


2.7 



218.600/228.600 
2.0 
0.8125 
a= 03 h 33 m 36?40 
<5=-36°08'25"7 
54/52 
100 x 50 
67, 28, 37, 4 



"Zenith opacity at 225 GHz measured at the Caltech Sub- 
millimeter Observatory next to the SMA. 
b Median system temperature toward NGC 1365. 
c Flux density of the calibrator. LSB/USB. 
d Integration time on the galaxy. 

c The full width at half maximum of the primary beam ap- 
proximated to be a circular Gaussian. LSB/USB. 

f Noise in the 12 CO, 13 CO, C ls O cubes and in the continuum 
data, respectively, used for Fig. [2] Those for line data are the 
ones in each channel. The velocity resolution is 10 km s -1 for 
12 CO and 20 km s -1 for other lines. 
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TABLE 3 

Measured and derived parameters 



parameter 


value 


unit 


5 ' 2 CO(2-l)( S S < 1 k P c ) 


4.3 x 10 3 


Jy km s 


s "co(2-i)( R s < 1 k P c ) 


3.8 x 10 2 


Jy km s _1 


°C l8 0(2-I)V iV g — 1 M'W 


1.1 x 10 2 


Jy km s 


max/i2 CO(2 _i) 


42 


Jy km s _1 arcsec~ 2 


D, 1£ LU(i-l) 


13 


K 


M mo i(R g < 1 kpc) 


9 x 10 s 


Mq 


max£ raol 


9 x 10 2 


Mq pc 2 


V sys (radio, LSR) 


1618 


km s -1 


R. A. dynamical center (J2000) 


03 33 36.35 


h, m, s 


Dec. dynamical center (J2000) 


-36 08 25.8 


o 1 II 


Td Y n(Rg = 1 kpc) 


3 x 10 7 


yr 


M iyn (R g < 1 kpc) 
M mol /M iyB (R g < 1 kpc) 


1 x 10 10 


M 


0.1 




mean 


8.7 




r.m.s. Rf 2 -« 


2.8 






1 x 10 s 




S m( ,i/S S FR(R g < 1 kpc) 


y 



NOTE. — Parameters related to gas mass are derived from 12 CO(2-l) 
data using A f (H 2 )// C 0(2-i) = 0.5 X 10 20 cm" 2 (K km s~' ) -1 and a factor 
of 1.36 for helium. See text for other assumptions and methods used to 
derive the parameters, and for their uncertainties. 



TABLE 4 

Association of CO Peaks with Known Sources 



sources association with a CO peak source parameters 



radio 


mid-IR 


optical 


In 


?b, 12 


hi 


?b, 13 


/is 


7b, 18 


S3. 5 cm 


S11.9 M m 


"'B 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Nuc. 


Ml 


Nucleus 


N 


N 


N 


N 


N 


N 


0.24 


510 ± 100 


17.00 


A 




SSC 10 


N 


N 


N 


N 


N 


N 


0.97 




21.42 




M2 


SSC 6 


N 


N 


N 


N 


N 


N 




< 20 


17.21 




M3 


SSC 3 


N 


N 


N 


N 


N 


N 




30±6 


17.20 


D 


M4 




Y 


Y 


Y 


Y 


N 


Y 


1.47 


40±8 




E 


M5 


b 


Y 


Y 


Y? 


Y 


N 


N 


0.80 


120 ±20 




F 






Y? 


Y? 


N 


N 


N 


N 


0.63 






G 


M6 




Y 


Y 


Y 


Y 


N 


Y? 


1.56 


140 ±20 






M7 




N 


N 


N 


N 


N 


N 








H 


M8? 




Y 


Y 


N 


Y 


N 


N 


0.58 






J 






N 


N 


N 


N 


N 


N 


0.37 







NOTE. — Col. (1) Radio sources iSandqvist et al. 1995; Forbes & Norris 1998). Col. (2) Mid- 
IR sources l Gall iano et all2005T) . Ml is the nucleus. Col. (3) Optical sources jKristen et alll997h . 
Cols. (4)-(9) Association with a peak in integrated intensity (/) and with a peak in brightness 
temperature (T b ). The subscripts 12, 13, and 18 indicat e l2 CO(2-l), "CO(2 -l), and C ,8 0(2-l), 
respectively. Col. (10) Flux density in mjy at 3.5 cm (Stevens et al. 1999). The numbers here 
differ from those in Table |S| because of the different methods used to estimate and subtract the 
background. Col. (11) Flux density in mjy in the '11.9 filter' of Galliano et al. 12005). Col. 
(12) B magnitude (Kri sten et alll997t) . 

'• Galliano^^lJ J2005) pointed out a conical feature near M4 in a HST image. See Fig. [5] 
fc Galli ano et alJ i2005) pointed out a double source near M5 in a HST image. See Fig. [3] 
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TABLE 5 

Properties of CO Hotspots with a Radio Counterpart 



R.A. 


Dec. 


velocity 


max7/ bI2 


hi 


AV U 


7b,12/7i,,13 




association 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


36?39 


18"7 


1560 


12.7 


42 


70 


7.3 ±0.6 


7.7 ±0.5 


D, M4 


36?61 


16'.' 


1540 


13.3 


34 


60 


7.8 ±0.7 


9.3 ±0.8 


E, M5 


36?77 


IS'.' 5 


1480 


10.2 


33 


80 


5.5 ±0.4 


8.1±0.7 


G, M6 


36?76 


23"6 


1510 


10.5 


32 


70 


9.9±1.4 


11.2±1.3 


H, M8? 


36?64 


28"0 


1590 


9.8 


37 


90 


12.9 ±2.2 


14.3 ±1.7 


F 



NOTE. — Cols. (l)-(2) Position. a=0 h 33 ra , 8 = -36°08' (J2000.0). Peak po- 
sitions measured in the peak temperature maps of 12 CO and 13 CO are averaged for 
the CO hotspots corresponding to D,E,G, and H. The radio position is used for F. 
Col. (3) Velocity at which the peak brightness temperature of 12 CO(2-l) emission 
is observed. The unit is km s~ . Columns from (3) to (6) are measured in the data 
cubes of 2."6 x l."8 x 10 km s _1 resolution. Col. (4) Peak brightness temperature of 
12 CO(2-l) emission in K. Col. (5) Integrated 12 CO(2-l) intensity in the units of Jy 
km s arcsec -2 . Col. (6) Full width at half maximum of the 12 CO(2-l) line. The unit 
is km s _1 . Col. (7) Ratio of 12 CO(2-l) to l3 CO(2-l) brightness temperature at the 
velocity where 12 CO(2-l) has the maximum intensity. The ratios in this and the next 
columns are measured in 3." 3 x 2"0 x 20 km s -1 resolution cubes made using visibili- 
ties in the same tt-v range. The error is 1 a. Col. (8) Ratio of 12 CO(2-l) to 13 CO(2-l) 
integrated brightness temperatures and its 1<t error. Col. (9) Spatial association with 
sources in other wave lengths. D through H are ra dio sources jSa ndqvist et all 1995: 
iForbes & Norrisll998 ) and M# are mid-IR sources iGalliano et aU2005l) . 



TABLE 6 

CO Hotspot- Cluster Associations 



NGC 1365 Antennae M83 (unit) 



distance" 


17.9 


21.5 


4.5 


Mpc 


log Lg_iooo M m d 


11.00 


10.84 


10.29 


L Q 


name(s) 


D/M4,E/M5,G/M6 


WS95-80/SGMC4-5 b 


M83-2 C 




S3.5cm 


1.8^4.2 d 


4.7 e 


4.8 f 


mjy 


^3.5 cm 


0.7-1.5 
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"The galaxy distance s and luminositie s are from ISanders et alJ 120031) . except that 
M83's distance is from T him et all 120 03 ) and its luminosity is scaled for the distance. 

b The names are from Whitmore & Schweizer 1 1995) and Wils on et alJ ^2000). The 
latter lists other names of the source. 



c In Table 2 of Turner & Hoi fT994T) . 

d These flux densities are measured from the 3.5 cm map of Morganti et al. 1 1999). The 
spatial resolution of the data is l"6x0."8. High-pass filtering is made for background 
subtraction using the AIPS task MWFLT with a 4"5 x 4"5 filtering box, to be consistent 
with the Antennae data. 

c The flux density of the source reported in lNeff & Ulvestad i200C|) . A l''4x l"0 reso- 
lution image was used after high-pass filtering . 

f Interporated from the 2 cm and 6 cm data in lTurner&HcliT994T) . The data resolution 
is 2"5 xl"6. No background subtraction is made and the source is assumed to be 
unresolved. 

sFlux in a 6" pixel on the source j Mirabel et al|1998tlWilson et aU20 00). 

h Flux in a 4" pixel on the source Tel esco et all!9 93). 

'Peak excess brightness temperature of 12 CO emission. The transition of CO is J=2- 
1 for NGC 1365 an d M83 and J=l-0 for the Antennae. The Antenna data are from 
IWilson et al. H2000l) and M83 data from Sakamot o'etai] 12004T) . 
^Integrated 12 CO intensity. 

k Peak surface density of molecular gas at the CO hotspot. It is estimated from the peak 
CO integrated intensity. Correction for inclination has not been made. 

'Resolution of the CO data, i.e., the FWHM of the observing beam at the distance of 
the galaxy. 
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FIG. 1 .— Spectra of 12 CO, I3 CO, and C 18 (J=2-l) emission in the central 25" of NGC 1365. For these spectra, the SMA data cubes also used for Fig. |3 
are first corrected for primary beam attenuation, then convolved to a resolution of 25" (FWHM), and sampled at the galactic center position of Sandavist et al. 
1 1982), which was used in many single-dish observations. The dashed line shows the systemic velocity of the galaxy, V sys (radio,LSR) = 1618 km s '. 
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Fig. 2. — CO line and continuum images of the center of NGC 1365. (a) CO(J=2-l) integrated intensity map. The n-th contour is at 9.1 xn 
Jy beam -1 km s~'. (b) Mean velocity map of the 12 CO(J=2-l) emission. Contours are in 30 km s~' steps, (c) 3 CO(J=2-l) integrated intensity 
map. Contours are at 2.7x[l,3,5,7,9, 11] Jy beam -1 km s~'. (d) Mean velocity map of the 13 CO(J=2-l) emission, (e) C' s O(J=2-l) integrated 
intensity map. Contours are at 3.6x [1, 3, 5] Jy beam -1 km s~' . (f) Continuum map at the effective frequency of 224 GHz, made by averaging 
the LSB and USB maps. Contours are at [-3,-2,2,3,4] x 4.0 mJy beam -1 (= la). Negative contours are dashed. The cross at the center of 
each image marks the active nucleus. The FWHM of the synthesized beam and a 500 pc scale are at the bottom left corner of each image. The 
dotted oval shows the mosaicked primary beam at its 50% level. The dashed line in the velocity maps is the line of nodes (P.A. = 220°). 
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FIG. 3. — 12 CO contours on an archival HST image taken with the F606W filter. Every 2nd contour in Fig. 0(a) is plotted on the logarithmically-scaled optical 
image. North is up, east is to the left. The HST image is shifted east by I'.'O from its header position to place the nucleus on its radio position. 




Radius [arcsec] Radius [arcsec] 

FIG. 4.— Radial distribution of 12 CO(2-l) emission around the galactic center measured in concentric circular annuli or circles in the galaxy plane, (left) The 
mean integrated intensities on the sky plane. The right axis is molecular gas surface density corrected for the inclination of the galaxy, (right) The cumulative 
flux of 12 CO(2-l) emission as a function of galactocentric radius. The right axis is the mass of molecular gas in units of 10 8 Mq . Both plots are from the dataset 
of 2"6 X 1"8 resolution. See $4. 1 I for the conversion of 12 CO data to gas surface density and gas mass. The two horizontal bars in each plot show FWHM of the 
elliptical beam along the galaxy's major axis (upper bar) and the minor axis (lower bar). The lengths of both bars are calculated in the galaxy plane. 
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FIG. 5. — Histogram of 12 CO(2-l) integrated intensity within 2.5 kpc of the nucleus in the galaxy plane. The vertical lines in the upper right mark the data for 
the CO hotspots in Tablel5l The moment map with 2'.'d X 1 " 8 resolution was used for this plot after corrected for the mosaicked primary beam. The area covered 
by the synthesized beam is 5 . 1 arcsec 2 on the sky and 5.1x1 0~ 2 kpc 2 on the galaxy plane. The top axis is the inclination-corrected surface density of molecular 
gas (S mo i) in units of M0pc~ 2 . See j|4.1l for the mass estimate from CO flux. The left axis is the surface area in the galaxy plane. Multiply the bottom axis by 
1.3 to obtain the total visual extinction in our line of sight in mag for N(H2)/Ay = 1 X 10 21 cm" 2 mag" 1 . 
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FIG. 6. — Position-velocity diagram of CO(2-l) emission along the line of nodes. The ra-th solid contour is at n X 134 mjy beam ; the lowest contour is 
at 2a. Negative contours are dashed and have the same absolute values as positive ones. Conversion to brightness temperature in K can be made by multiplying 
by 5.2. The diagram is corrected for the mosaicked primary beam. The abscissa is measured from the active nucleus. The horizontal dotted line indicates the 
systemic velocity of 1618 km s . The black rectangle at the bottom right corner represents a resolution element. 
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FIG. 7. — (left) Comparison of 12 CO(2-l) and 13 CO(2-l) line intensities in the center of NGC 1365. The line intensities are sampled from data cubes of 
matching spatial and velocity resolutions, 3"3 X 2"0 and 20 kms~'. The sampling spacing is lf'2. Only those points where both lines are detected above 
3cr are plotted. The la noise is 0.12 K and 0.09 K for 12 CO(2-l) and 13 CO(2-l), respectively, at the center of the mosaics. The noise varies across the maps 
corrected for primary beam, and so does the cutoff threshold in kelvin. (right) Histogram of 12 CO(2-l) to l3 CO(2-l) intensity ratio. The one in black is from 
the entire data. The blue and red ones are from the near and far side, respectively, and within the galactocentric radius of 1 kpc (= 12"). Data from outside the 
circumnuclear disk are plotted in green. 




I 3.3" x 2.0" 500 pc 



t i i i i i i i i i i i i i i i i i i i J 

20 10 -10 -20 

R.A. offset [arcsec] 

FIG. 8. — Spatial distribution of 12 CO(2-l) to 13 CO(2-l) brightness temperature ratio. Ratios of different velocities are averaged at each position. The + sign 
marks the active nucleus and the dashed line is the line of nodes; the near side is to the northwest. The dotted ellipse is at the galactocentric radius of 1 kpc (= 
12"), and the boxes and crosses mark the radio and mid-IR sources, respectively, in Fig.l9l 
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FIG. 9. — A guide map of radi o, mid-IR, and optical sources in the center of NGC 136 5. Sources A, D, E, F, G, H, and J are rad io hotspots. Their positions are 
shown as diamonds i Sandqvist et al. 1995, A„|, s = 2. 6 cm), squares I Stevens et al. 1999, A„|, s = 3.5,6 cm), and circles i Thean et al. 2000, A (, s = 3.5 cm). Sources 
M2 through M8 sh own as X sym bols are mid-IR sources of Galliano et al. (2005). Super star clusters (SSCs) 3, 6, and 10, shown with plus signs, are optically 
identified by Krist en et alj 119971) . SSC 3 and SSC 6 are the two brightest (Mb = -16.6 mag) clusters in the region. Positions are measured with respect to the 
Seyfert nucleus; see Table^for its absolute coordinates. The background is our map of I2 CO(2-l) integrated intensity to guide the eyes. 
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FIG. 10. — Radio and mid-IR sources plotted on the maps of CO integrated intensity (left) and peak brightness temperature (right). Squares are seven 
radio peaks observed in centimeter wavelengths i Stevens et al. 1999; Sandqvist et al. 1995; Thean et al. 2000) and crosses are mid-IR sources observed by 
IGalliano et al. 1 2005). The background images are integrated intensities in Jy beam -1 km s -1 for the left column and brightness temperatures in K for the right. 
Each map is corrected for the primary beam attenuation. The + sign at the center of each image marks the galaxy nucleus. The synthesized beam and a 500 pc 
scale are at the bottom left corner of each panel. 
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FIG. 1 1 . — Spectra of CO hotspots. At the top of each panel are the radio name of the hotspot and the line name. The spectra are sampled from the data cubes 
used for Fig.l2lat the positions in TablelSlafter primary-beam correction. Intensities are in brightness temperatures. The beam size, and hence the degree of beam 
dilution, are different for the three lines. 
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FIG. 12. — (a) Spectra of 12 CO near an apocenter of the circumnuclear ring. The double-peak profiles indicate the motion of gas roughly along the circumnuclear 
ring, (b) Sampling positions of the line profiles, shown as a group of small plus signs. The background image is the 12 CO integrated intensity, and the contours 
arc the 12 CO mean velocity in 20 km s steps, (c) Illustration of the gas motion. The double-peaked profiles in the sampling region (= the gray square) are due 
to the distinctively different velocities in our line of sight between the gas approaching on the ring to the apocenter and the gas leaving it along the ring. The gas 
streamline on the ring can be oval similar to the dotted ellipse, or it may have a cusp at a shock. In either case, its major axis makes an oblique angle with the 
stellar bar (< 90° measured in the direction of the galaxy rotation, or clockwise), indicating that the gas streamline is .t2-like. Around the apocenter, the gas flow 
on the molecular ring and another flow roughly along a leading-edge gas lane on the bar (almost) converge, likely leading to frequent cloud collisions. The major 
and minor axes of the galaxy are shown as the dashed lines and the stellar bar as the dot-dashed line. 
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FIG. 13. — NGC 1365 in the Schmidt-law plot of Kennicutt 1 1998, Fig. 5) for starburst nuclei. The gas surface density and the star formation rate are averaged 
over the circumnuclear starburst region of each galaxy. The size of the starburst region is individually determined and is different for each galaxy. It is 2 kpc in 
diameter for NGC 1365 and is taken from K98 for others. The gas surface densities are normalized with xqo= 2.8 only to make a fair comparison. Asterisks 
are mergers and diamonds are non-merging galaxies. The plu s symbo ls are for two distant (> 100 Mpc) galaxies that are likely mergers judging from their large 
far-IR luminosities (log(Lg_iooo;im/£Q) > 11.80; Sanders et al. 2003) but lack morphological evidence for merging. 



